Qubit Disentanglement and Decoherence via Dephasing 
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We consider whether quantum coherence in the form of mutual entanglement between a pair of 
qubits is susceptible to decay that may be more rapid than the decay of the coherence of either 
qubit individually. An instance of potential importance for solid state quantum computing arises 
if embedded qubits (spins, quantum dots, Cooper pair boxes, etc.) are exposed to global and local 
noise at the same time. Here we allow separate phase-noisy channels to affect local and non-local 
measures of system coherence. We find that the time for decay of the qubit entanglement can be 
significantly shorter than the time for local dephasing of the individual qubits. 

PACS numbers: 03.65.Ta, 03.65.Yz, 03. 67. -a 



I. INTRODUCTION 



Quantum coherence is fragile in the sense that when a 
quantum system is in contact with an environment with 
many degrees of freedom, a quantum superposition of 
a pointer basis inevitably undergoes a decoherence pro- 
cess pi Q , which is characterized by the decay of the off- 
diagonal elements of the density matrix in the given basis. 
The fragility of quantum coherence is responsible for the 
lack of observation of Schrodinger Cat states in everyday 
life 0- Quantum entangled states, as special forms of 
coherent superpositions of multi-partite quantum states, 
have recently been recognized as a valuable resource that 
is of crucial importance in realization of many quantum 
information protocols, such as quantum cryptography 
El Q- quantum teleportation [ 1, s u per dense coding 
J8j, and quantum computation [(J, flCX fill Appli- 
cations of interest, such as issues in quantum informa- 
tion processing (QIP) 0,0], have triggered extensive 
research aimed at controllin g qu antum disentanglement 
CI 111 III 13 13 13 El HM E|. Apart from the im- 
portant link to QIP realizations, a deeper understanding 
of entanglement decoherence is also expected to lead to 
new insights into the foundations of quantum mechanics 

mm 

The distinctions between "local" and "non-local" de- 
coherence can be important in any application in which 
coupled qubits may be remote from one another but nev- 
ertheless are coupled or jointly controlled through com- 
mon external influences such as electromagnetic fields 
(microwave, optical, magnetic, etc.). In this paper, we 
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focus on these distinctions in the simplest case, a pair 
of qubits (e.g., spins, quantum dots, Cooper pair boxes, 
etc., which we will refer to as spins) embedded in in a 
solid state matrix and subjected to the most elemen- 
tary relaxation mechanism, pure dephasing. However, 
we include more than one source of dephasing noise and 
allow the different noise sources to produce dephasings 
with different time scales, in order to mimic a relatively 
general practical situation. With these simple elements 
we address a few general questions: How does dephas- 
ing affect entanglement as well as local coherence? More 
concretely, how rapidly does disentanglement occur com- 
pared to the off-diagonal decay rates of a single qubit 
density matrix? and how is the disentanglement rate re- 
lated to those other coherence decay rates? 

Our results, obtained relatively directly through ex- 
plicit derivation of the 12 relevant Kraus operators for 
the several decoherence channels, are as follows: We show 
that an environment that causes pure dephasing can af- 
fect entanglement and coherence in very different ways. 
As a consequence, we show that disentanglement gen- 
erally occurs faster than the decay of the off-diagonal 
dynamics of either the composite two-qubit system or 
an individual qubit. One manifestation of the difference 
is that decoherence of a composite two-qubit state can 
be incomplete, but entanglement will be completely de- 
stroyed after a characteristic time. Our findings, while 
not entirely in contradiction to one's general intuition, 
may still come as a surprise. 

Specifically, we will investigate non-local and local as- 
pects of coherence decay when two qubits are subject 
to different dephasing environments. This situation may 
naturally arise in a long distance quantum communica- 
tion, local manipulation of qubits, or in a local quantum 
measurement [3 13 |3 13 With application of the 
quantum map approach to our two-qubit model, we ex- 
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plicitly show that entanglement, measured by Wootter's 
concurrence [3(J, generally decays faster. We show that, 
by altering environmental parameters, the time for dis- 
entanglement can be significantly decreased compared to 
the usual time for complete coherence decay. The present 
work extends an earlier result where two qubits are as- 
sumed to interact collectively with a single reservoir |3l| . 

The paper is organized as follows. We introduce a two- 
qubit dephasing model and describe it in the language of 
noisy quantum channels in Sec. [H] and Sec. 11111 respec- 
tively. The explicit solutions to the model are given in 
Sec. IIVI In Sec. |3 we study the dephasing processes for 
the channel that acts on both qubits and for the channel 
that acts on only one qubit. We investigate disentan- 
glement in Sec. EH In Sec. EH1 

we discuss coherence 
decay in terms of transmission fidelity. We comment and 
conclude in Sec. IVIIII 




II. TWO-QUBIT SYSTEM WITH THREE 
RESERVOIRS 



FIG. 1: The model consists of two qubits A and B collectively 
interacting with the stochastic magnetic field B(t) and separately 
interacting with the stochastic magnetic fields b^(t) and bg(t). 



We consider two qubits A and B that are coupled to 
a noisy environment both singly and collectively. One 
specific realization would be a pair of spin-i particles 
embedded in a solid-state matrix and subject to random 
Zeeman splittings, for example from random extrinsic 
magnetic fields. Analogously, one may envision a pair of 
polarized photons travelling along partially overlapping 
fiber-optic links which are gradually de-polarized due to 
random birefringence. In the first case the relaxation de- 
velops in time and in the second case in space. The bare 
essentials of the physics can be described by the following 
interaction Hamiltonian (which takes h — 1 and adopts 
spin notation): 



H(t) = -^(B(t){v?+*?) + b A (t)ir* + b B (t)iTf), (1) 

where p is the gyromagnetic ratio and Z?(f), 6s(i) 
are stochastic environmental fluctuations in the imposed 
Zeeman magnetic fields. Here crf ,B are the Pauli matri- 
ces: 



A.B 



1 

-1 



(2) 



and we take as our standard 2-qubit basis: 



\1)ab = 
\3)ab = 



)ab, \2)ab = 
)ab, |4)ab = 



)AB, 

)ab, 



(3) 



where | ± ±)ab = <H) |±)b denote the eigenstates of 
the product Pauli spin operator o^®o® with eigenvalues 
±1. 

For simplicity, we assume that the stochastic fields 
B(t),bA(t),bB(t) are classical and can be characterized 
as statistically independent Markov processes satisfying 



< B(t) > 

< B{t)B{t') > 

< bi{t) > 

< bi(t)bi(t') > 



0, (4) 

j-S(t-t'), (5) 

0, (6) 

^S(t-t'), i = A,B (7) 



where < ... > stands for ensemble average. Here T is the 
damping rate due to the collective interaction with B(t), 
and Ti are the damping rates of qubit i (i = A, B) due 
to the coupling to the fluctuating magnetic field bi(t). 
The white-noise properties in JSJ and J7J ensure that the 
two-qubit system will undergo a Markov evolution. More 
realistic and more general models are easy to imagine, 
but for the sake of simplicity we will restrict ourselves to 
this simple situation. 

Hamiltonian represents a class of models of open 
quantum systems that describe a pure dephasing pro- 
cess, which can arise from interaction with Bosonic en- 
vironments, independent of whether they are associated 
with quantized thermal radiation fields or phonon fields 
in solids or other specific physical contexts. 

The dynamics under the Hamiltonian Q can be de- 
scribed by a master equation, but for our purpose we find 
it more convenient to use the Kraus operator-sum repre- 
sentation The advantage of the operator-sum rep- 
resentation will become clear in Sec. lVIl The density ma- 
trix for the two qubits can be obtained by taking ensem- 
ble averages over the three noise fields B(t), b A it), bsit): 



p{t) =«< p s tit) >>>, 



(8) 
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where the statistical density operator p s t(t) is given by 



Pat (t) = U(t)p(p)U*(t), 
with the unitary operator 



U(t) = exp 



-i [ dt'H{t') 
Jo 



(9) 



(10) 



III. 



KRAUS REPRESENTATION FOR NOISY 
CHANNELS 



We can describe a decoherence process in the language 
of quantum channels |3^ |. We take a quantum channel 
to be a completely positive linear map that acts on the 
quantum state space of a system of interest. Let £ be a 
quantum channel that maps the input state p ln = p(0) 
into the output state p ut = p(t). It is known that the 
action of £ can be characterized by a (not unique) set 
of operators called Kraus operators [Til Il2l |3f| . For 
any initial state, the action of quantum map £ is given 
by 



N 



p(t) = £{p) 



£4(*)p(o)i^(i), 



where are Kraus operators satisfying 



(ii) 



(12) 



for all t. Eq. (|llfl is often termed the Kraus (or operator- 
sum) representation in the literature (for example, see 
[Til [l2|p. The operators contain all the information 
about the system's dynamics. Unitary evolution is a spe- 
cial case with only one Kraus operator; otherwise the 
channel describes a non-unitary process associated with 
damping and decoherence. It can be easily seen that any 
Kfj, projects pure states into pure states, but the collec- 
tive action of different K^s will typically transform a 
pure state into a mixed state. 

The most general solution (JSJ can be expressed in 
terms of twelve Kraus operators (under the assumption 
that the initial density matrix is not correlated with any 
of the three environments): 

2 3 

p(t) = £(p(0)) -EE D l E j F iP(°) F i E i D k, (13) 

i,j=l k—1 

where the Kraus operators describing the interaction 
with the local environmental magnetic fields bA{t), &s(i) 
are given by 



E 1 = 



Fi=It 



1 

7a(<) 
1 
7fl(*) 



I, E 2 = 



,F 2 =I 





u A (t) 



w B (t) 



I, (14) 
(15) 



and the Kraus operators that describe the collective in- 
teraction with the environmental magnetic field B(t) are 
given by 



Do = 



Da 



hit) 
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(16) 



(17) 



(18) 



\0 uj 3 {t)j 



Note that the parameters appearing in (|L4|) — (|LH|) are 
given by 



lA(t) 
UJ A (t) 

7(i) 
wa(t) 

wa(t) 



= p -t/2T, 



2 , 7s(*) 



(19) 

\T\ -e-t/Tf, uj B (t) = Vl-e-t/^, (20) 

-e-*/ T 'y/T^, 



-t/T 2 



^(l_ e -t/T 2 )( 1 _ e -2t/T 2 ) j 



(21) 
(22) 
(23) 



where T 2 = 1/r is the phase relaxation time due to 
the collective interaction with B{t) and T 2 = 1/Ta and 
T 2 B = 1/Yb are the phase relaxation times for qubit A 
and qubit B due to the interaction with their own en- 
vironments bA{t),bB{t), respectively. We will no longer 
write the time-dependent arguments of the 7's and w's 
as a notational simplification. 



IV. EXPLICIT SOLUTIONS, VARIOUS 
DENSITY MATRICES 

Three cases can be distinguished: (1) The case that the 
two qubits only interact with their local environments, 
namely B(t) = 0; and (2) the case in which only one 
qubit is affected by a local environment, that is, B(t) = 
and one of 6>i(f), 6_b(£) is also zero. We call the above two 
cases the two-qubits local dephasing channel and the one- 
qubit local dephasing channel, respectively. (3) The case 
that the two qubits only collectively interact with B(t). 
We call this case the collective dephasing channel. We 
will present these three cases separately in what follows. 



A. Two-qubit local dephasing channel Sab 

When the pair of qubits only interact with the local 
magnetic fields, i.e. B(t) = 0, then the phase-noisy chan- 
nel, denoted by £ab, is characterized by the following 
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four composite Kraus operators: 



C. Collective dephasing channel £d 



Mi = EiF u M 2 = EiF 2 , 
M 3 = E 2 Fx, M 4 = E 2 F 2 , 



(24) 
(25) 



where the fundamental Kraus operators E% , E 2 , Fx , F 2 
are defined in (|14p. and (|15|) . Therefore, the effect of 
the quantum channel Sab on an initial state p can be 
described as follows: 

4 

P (t) = £m(p(0)) = X! M >(0)^- (26) 
The explicit solution in the standard basis is given by 



p(t) = £ab(p(0)) 

( Pll IBPII 

IB Pll P22 

7AP31 7A1BP32 

\lAlBPil 1AP4 2 



1AP13 7A7SPl4\ 

1A1BP 2 3 1AP2A I 

P33 7-BP34 I 

7-BP43 P44 / 



(27) 



where and jb are defined in (|19|l and the notation pij 
stands for the initial state pij(0) 

One immediate conclusion that can be drawn from 
Eq. (|2*?jl is that the two-qubit dephasing channel has an 
effect on all the matrix elements except for the diagonal 
elements. Thus, decoherence takes place in every super- 
posed state of the basis elements (0. Furthermore, it 
will turn out later in Sec. IVfl that this channel generally 
results in faster disentanglement. 



B. One-qubit local dephasing channels £a and Sb 

A simple case arises if the channel is assumed to affect 
only one qubit, say qubit A. In this situation, we call the 
channel a one-qubit local dephasing channel, and denote 
it by Ea- The one-qubit dephasing channel Ea can be 
described by the Kraus operators Ei,E 2 . Explicitly, for 
any initial state p(0), the output state under this channel 
is given by 



P (t) = Ea(p(0)) = Y, E Ip(°) E » 



( Pn 

P21 



Pll 
P 22 



7A/913 7A/9l4 N 
7AP23 7AP24 



1AP31 7AP32 
\lAPil 7A/042 



P33 
P43 



P34 
P44 



(28) 



Obviously, one can equally well define in terms of Kraus 
operators F\ , F 2 the one-qubit local dephasing channel 
£b, that has effect only on qubit B. 

Since the one-qubit dephasing channels only affect a 
single qubit, the coherence of the other qubit is not af- 
fected during processing. Consequently, the coherence 
of the composite two-qubit system cannot be destroyed 
completely. The details of the loss of coherence and 
entanglement under the channels defined above will be 
given in Sec. El and Sec. I VII 



If the local stochastic fields fe^i), frs(i) are negligible 
or switched off, we obtain a model that describes the 
two qubits collectively interacting with a common envi- 
ronment 0,H2- In this case, the solution p(t) can be 
cast into the operator-sum representation in terms of the 
three Kraus operators D^: 



3 

P (t) = £ D (p(0)) = ]T Dl(t)p(0)D^t), (29) 



where D^, (p = I, 2, 3) are given in JTHJ), (JT7|l and (|T%|) . 

The explicit solution of (|29|l in the standard basis © 
can then be expressed as 



p(t) = £ D (p(0)) 



Pll 1P12 1P13 7 Pl4 N 

7P21 P22 P13 7P24 

7P31 P32 P33 7P34 

V7 4 /°41 7P42 7P43 Pa 



(30) 



where 7 = e~*/ 2T2 is defined in 12I|) . 

As seen from (|30|) , the collective dephasing channel af- 
fects the off-diagonal elements P14, pn more severely than 
it affects the other off-diagonal elements. As a result, we 
have already shown in 31] that, for a set of entangled 
pure states containing P14, /341, disentanglement proceeds 
with rates that are faster than the dephasing rates for 
qubit A or qubit B. In the case of collective coupling, the 
stochastic magnetic field B{t) always affects the phases 
of qubit A and qubit B in the same way, so one expects 
that for some states the random phases may cancel each 
other out. Here such a cancellation due to the under- 
lying symmetry of the interaction Hamiltonian can be 
expected and has been demonstrated in various contexts 
[33. For the situation considered here, the cancellation 
leads to the existence of the disentanglement-free sub- 
spaces spanned by a basis consisting of \2)ab, |3)ab- 
The entangled states belonging to these subspaces are 
robust entangled states |3lj . 



V. LOCAL DEPHASING AND MIXED 
DEPHASING 

First, we study the coherence decay of a single qubit 
under the two-qubit dephasing channel £ab defined in 
(|2*7Jl . The local dephasing rate of the qubit can be de- 
termined from the reduced density matrices s A and s B , 
obtained from the density matrix (|27|l in the usual way, 
that is, 



Tr B {p} and s B = Tr A {p}- 



(31) 



The reduced density matrices for qubit A and qubit B 
are thus obtained as: 



s A (t) 



' 'pn(t) + p 22 {t) Pl 3 {t) + p 2i (t) s 

K P31 (t) + P42 (t) P33 (t) + P44 (t) , 



(32) 
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and 



'(*) = 



Pll(t)+p33(t) Pll{t) + P3i(t) 
P2l (t) + P43 (t) P22 (t)+p44(t) 



(33) 



The dephasing rates are determined by the off-diagonal 
elements: 



4« 



pl 3 {t) + P2i(t) =~/As£ 2 (0), 



and 



=Pi2(*) + P34(i)=7Bsf 2 (0). 



(34) 



(35) 



Thus, the dephasing times denoted by ta,b for qubits A 
and B, respectively, can be read off from the expressions 
for 7^ and jb- That is 



ta 



2 



^— , and tb = ^— 



2 



(36) 



Similar to local dephasing processes, the decoherence 
rate for the composite two-qubit system is characterized 
by the decay rates of the off-diagonal elements of the 
density matrix (|27|l . As seen from l|27|l . the two-qubit 
off-diagonal elements Pij(i < j), are given by 



pij(t) 



-Tat 



so the coherence pij (t) decays on a timescale 



7"de 



r, 



(37) 



(38) 



where r i2 = T 34 = 1^/2, Ti 3 = T 24 = T A /2,T 1A = 

r 23 = (r A + r B )/2. 

We denote r as a time scale on which all the off- 
diagonal elements of the density matrix i(2"7|) disappear. 
The corresponding decoherence rate 1/r is called here 
the mixed dephasing rate in order to distinguish it from 
the local dephasing rates of a single qubit. It is evident 
that the mixed dephasing rate 1/r is determined by the 
slower decaying elements [3(J, and in general, we have, 



T > T A ,T B - 



Namely, and this is a key point, the mixed dephasing 
rate 1/r is not shorter than the local dephasing rates for 
qubit A or qubit B defined in i|36|) . As will be seen in the 
next section, this is in contrast with the disentanglement 
rate, which can be faster than the local rates 1/ta and 

1/tb. 

Finally, we study dephasing processes under the local 
one-qubit dephasing channels Ea and £b- Since the local 
dephasing channels only affect one qubit and leave the 
other intact, hence some coherence may always exist in 
the composite two-qubit state. Let us consider the chan- 
nel £a and a superposed state of just three members of 
the standard basis, 



MB 



U\1)ab + \2)ab + H)ab)- (40) 



y/3 



Obviously, from (|28|l . we see there is no way that the 
channel £a can destroy the off-diagonal elements pi2, p 2 i, 
so £a cannot completely destroy the coherence of the 
composite two-qubit state l|40f) . Next, we look at the 
dephasing processes for an individual qubit. Obviously, 
the coherence of qubit B is not affected by £a- That is, 



'(*) = AO)- 



(41) 



Notice that s A (t) under the channel £a is given by H32fl . 
Hence, as expected, the coherence of qubit A will be de- 
stroyed by the channel £a on the dephasing time scale 
ta- A similar analysis applies to £b- 



VI. ENTANGLEMENT DECAY 

To describe the temporal evolution of quantum en- 
tanglement we have to measure the amount of entangle- 
ment contained in a quantum state. Since the entangle- 
ment decoherence processes are typically associated with 
mixed states, we will use Wootter's concurrence to quan- 
tify the degree of entanglement [3(j ■ The concurrence 
varies from C = for a disentangled state to C = 1 for a 
maximally entangled state. For two qubits, the concur- 
rence may be calculated explicitly [20] from the density 



matrix p for qubits A and B: 



C(P) 



(o, 



(42) 



where the quantities Xi are the eigenvalues of the matrix 



(43) 



arranged in decreasing order. Here p* denotes the com- 
plex conjugation of p in the standard basis © and a y is 
the Pauli matrix expressed in the same basis as: 



A.B 



-i 

1 



(44) 



(39) For a pure state the concurrence (|42|l is reduced to 



(45) 



A. Entanglement decay under two-qubit dephasing 
channel £ab 



How does the concurrence behave under a noisy chan- 
nel? From l)45p. for any entangled pure state 

I*) = oi|1)ab + a 2 \2) AB + a 3 \3) AB + a 4 \4:) AB) (46) 

where X^=i l a «| 2 = 1j the concurrence of the pure state 



ffify is given by 



C(|*)) = 2\ ai a 4 - a 2 a 3 \. 



(47) 
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It follows from l|27(l that the concurrence (|47|l will always 
decay to zero on a time scale determined by the dephasing 
time t defined in H39JI . In fact, let us first note that the 
density matrix of l|46l) will approach a diagonal matrix 
at the two-qubit dephasing rate 1/r. Since a diagonal 
density matrix only describes classical probability, then 
the concurrence must be zero. 

Moreover, we show in what follows that all entangled 
(possibly mixed) states decay at rates that are faster than 
the dephasing rates of an individual qubit. For this pur- 
pose, we first note that the concurrence C(p) is a convex 
function of p |3Cj : that is, for any positive numbers and 
density matrices p^, ([i = 1, ..,n) such that = 1> 

one has 



fj,—i fi—i 



(48) 



Applying the above inequality to l|27|) . we immediately 
get 

4 

C(p(t)) < ]T C(Mtp(0)M M ), (49) 

where M M are defined in (|24Jl and H25[l . Now we consider 
a typical term C{Mj 1 p{0)M l _ l ) in (gUl and denote it by 

p out = Mtp(0)M M . (50) 
Substituting (|5H|l into l|4^|l . we have, 

Q = Pout(^®<Kut(<<H>0 (51) 
= MlpM^®^)Mlp*M^®a^). (52) 

Notice that g has the same eigenvalues as g' 

g' = pM^a* ® a*)MTp*M;(a£ ® <f )Aft. (53) 

Also note that the Kraus operators defined in (|14|l satisfy 
Mj = — M* . Hence it is easy to check that 

M,K ® af)M, = { & (54) 

Then Eq. l(53l becomes: 

J> _ J iWbpK ® <K« ® O if M = 1, 







if 



^i: (55) 



By using the definition of concurrence (|43|) and the in- 
equality l|l9")l. we get: 



W)) < 7 a7bC(p(0)). 



(56) 



Hence, the entanglement decay time r e can be identified 
as 



11 1 

T e TA TB 



(57) 



where the local dephasing times ta, t~b are given in l|36|l . 
In this way we immediately see that the entanglement 
decay times for the two-qubit channel are shorter than 
the local dephasing times ta,tb and hence are shorter 
than the mixed dephasing time r as well. 

By suitably choosing the channel parameters Ta,^b, 
one can achieve 



r e < r, 



(58) 



where r is the mixed dephasing time (|39[) . We have 
shown here that the decay time for the entanglement of 
the two-qubit system can be much shorter than the mixed 
dephasing time r. An interesting case occurs when the 
phasing damping rates for qubit A and qubit B are as- 
sumed the same: ta — tb, then from (|57|l . we have 



Li 

2 



Tb 
2 



(59) 



This relation reminds us of the well-known relation be- 
tween the phase coherence relaxation rate T 2 and the di- 
agonal element decay rate T\ in open quantum systems 

To see how the phase-noisy channels affect entangle- 
ment in a more explicit way, it would be instructive to 
find some examples in which the concurrence can explic- 
itly be calculated. For this purpose, let us consider, for 
example, the following two classes of almost arbitrary 
bipartite pure states described by: 

|<£i) = ai|l) A s +a 2 \2) AB +0414)^, (60) 
1^2) = ai\l) A B + a 3 \3) A B + a^tyAB, (61) 



and 



1^2) 



oi|1)ab 
02|2)ab 



■ a 2 \2)AB 
' «3 1 3) AB 



■ C13\3}aB, 
- a,i\A) A B- 



(62) 
(63) 



It will turn out in what follows that the concurrence for 
those two classes can be calculated exactly. 

For the pure state l|6T?|) . the concurrence is C(|</>i)) = 
2|ai<24|. Then the density matrix with the initial entan- 
gled state jlSTj|) at t is given by 



Pit) 



£ab(p) 

( M 2 

lBa 2 al 




7 B aia2 7A7s a i a 4 N 
|a 2 | 2 7,40204 
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(64) 



KlAlBa^al ^ A aia* 2 |a 4 | 

The non-zero eigenvalues of the matrix g defined in O 
are 

Ai = (l+7 J 47_ B ) 2 |aia 4 | 2 , A 2 = (l-7A7s) 2 |aia 4 | 2 . (65) 

Then the concurrence at t is given by 

C(p(t)) =2 7A 7b|oi04|. (66) 

Similarly, as another illustration, we may calculate the 
entanglement degree of the entangled states (|62|l . The 
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density matrix with the initial entangled state l|62[l at t 
is given by 



p(t) = £ab{p) 



( kl 2 

lBa 2 a\ 

V o o 



7saia 2 
|a 2 | 2 



7Aaia^ 0\ 

lAlBa>2al 

|a 3 | 2 

0/ 



The concurrence at t can be obtained 



(67) 



(68) 



From the above two examples, we see that the concur- 
rence is determined by the fast decaying off-diagonal ele- 
ments P14 and p 2 3- The examples also demonstrate that 
the bound in (|56|l is the minimal upper bound. 



B. Entanglement decay under one-qubit dephasing 
channels £a or £b 

Perhaps the best way of seeing the difference between 
the disentanglement and decoherence is to look at the 
effect of the one-qubit channel on a two-qubit system. 
Similar to the two-qubit dephasing channels, it can be 
shown that under the one-qubit local dephasing channels, 
the concurrence is determined as follows: 



C(£a(p(0)) < 7A<7(P(0)), 
C(£ b (p(0)) < 7flC(p(0)). 



(69) 
(70) 



Thus, we see that the one-qubit local dephasing channels 
can completely destroy the quantum entanglement after 
the dephasing times ta or tb ■ To be more specific, let 
us consider the following entangled pure state: 



1 



J) 



AB 



+ \3)ab + |4) 



AB 



(71) 



If we apply the local dephasing channel £ a on p = \<f>){<f>\ 
then the output state is given by 



p(t) = £a(p) 



1 n 1A ia 

3 u 3 3 



1 1 

in 1 ? 

3 3 



(72) 



The concurrence of (|72() can be obtained 



COM) = ~e~ t/2Ti 



(73) 



This means that disentanglement will proceed exponen- 
tially in the local dephasing time ta- As seen from Sec.lVl 
however, the coherence of the composite two-qubit state 
\<j)) cannot be destroyed completely. In general, one-qubit 
local dephasing channels £a and £3 cannot destroy the 



coherence between different elements of the basis l|46ll . 
That is 



4 

p(*) + E 

i=l 



\i){i\, as t 



(74) 



Namely, the decoherence process for the composite two- 
qubit system becomes frozen after the local dephasing 
times ta or tb- 

In summary we have shown that the effect of the de- 
phasing channels on entanglement and quantum coher- 
ence may differ substantially. In general, the entangle- 
ment of two qubits decays faster than the quantum co- 
herence. In particular, for the channels £a and £b, we 
have shown that for some entangled states the local de- 
phasing channels cannot completely destroy the coher- 
ence between the qubits. However, entanglement as a 
global property can be completely destroyed by the one- 
qubit local dephasing channels. 



VII. INPUT-OUTPUT FIDELITY AND 
CONCURRENCE 

It is interesting to look at the coherence and entan- 
glement decays by comparing input-output fidelity and 
concurrence. The fidelity is meant to measure the "differ- 
ence" between an input state and an output state under 
a noisy quantum channel. As such, it has been used to 
measure the loss of coherence of a quantum state when 
the state is subject to the influence of an environment. 
The fidelity can be defined as 



F(pi n ,Pout) 



(75) 



where the maximum is taken over all purifications 
IV^) of p m and p out [3^, |3|j. For a pure initial state 
Pin = IV'KV'Ij ^ can simply be written as 



F(pi n , Pout) 



lA>ut#)- 



(76) 



In the case of the one-qubit dephasing channel £a 
defined in Q28J1. for the general entangled states l|46|l . 
Eq. I)76|) can be written as 



f(p,£a(p)) 



4 

E 

1=1 



27^ (|aia 4 | 2 
27A (|aia 3 | 2 
2(1 



aia 2 



|a 2 a 3 | 
\a 2 a 4 \ 

|2\ 



|a 3 a 4 | 



(77) 



For example, we consider the entangled pure state l)46|l 
with coefficients 



01 = «2 = a-3 



1 1 

2' ai = —2- 

Then it is easy to check that, for all times t, one has 

1 



(78) 



F{ P ,£ A {P)) > 



(79) 
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Hence we again see that some amount of coherence is 
always preserved. However, as shown in the previous 
discussions, it is obvious that the disentanglement will 
be complete; that is, C decays to zero on the scale of the 
local dephasing time ta- 

C( P ) = 1, C{E A { P )) = 0. (80) 

VIII. CONCLUDING REMARKS 

In this last section we provide some comments on dis- 
entanglement, mixed dephasing processes, and local de- 
phasing processes for an individual qubit. 

It is evident from sections Ivl and IvTI that in general en- 
tanglement is more fragile than local quantum coherence. 
In the case of the examples H66|) and B68JI . we see that 
quantum entanglement is determined by the off-diagonal 
elements pu, P23 and their conjugates. In our two-qubit 
model, these elements are affected more severely than 
other off-diagonal elements. So the disentanglement nor- 
mally proceeds in a rate which is faster than the dephas- 
ing rate, the latter being generally determined by the 
slow-changing off-diagonal elements pi2,Pi3 etc. 

We remark that our simple approach to decoherence in 
terms of phase damping channels can easily be extended 
to other decoherence channels such as amplitude damp- 
ing channels. Although the disentanglement rates are 
dependent on the choice of the measure of entanglement, 
it seems that the fast disentanglement demonstrated in 
the paper is of generic character. 

Since the local manipulations of a single qubit can- 
not increase the entanglement between two qubits, so 
the quantum map 1|27[1 cannot create any entanglement. 
This is in contrast to the collective dephasing channel 
£r>- It is known that entanglement may be created due 
to the interaction with a common heat bath [3^ . 

In connection to the foundations of quantum theory, a 
good understanding of entanglement decoherence might 
provide a more detailed picture of the transition from 



quantum to classical (or semiclassical) dynamics. As 
suggested in this paper, a quantum system with several 
constituent particles, if initially in an entangled state, 
may undergo several stages of transition, including the 
"collapse" of entanglement, and then the "collapse" of 
coherence among individual particles. This is a topic to 
be addressed in future publications. 

Summarizing, our primary interest in this paper is 
to investigate the fragility of quantum entanglement. 
Through simple two-qubit systems we have shown that 
the most commonly assumed (dephasing) model of noisy 
environments may affect entanglement and quantum co- 
herence in a very different manner. While preserved in 
the case of collective interaction, the permutation sym- 
metry of the two qubits is violated by local environ- 
ment interactions. As a result, there are no entangle- 
ment decoherence-free subspaces for the local dephasing 
models, which means that entanglement for all entangled 
states will decay into zero after some time. Our results 
suggest that, when spins are coupled to inhomogeneous 
external magnetic fields, disentanglement and dephasing 
may proceed in rather different time scales. It means 
that, even in the situation that quantum coherence has 
relatively long relaxation times, entanglement may be 
rapidly suppressed. This issue is certainly of importance 
in realizing qubits as spins of electrons or nuclei in solids. 
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